Kinetic Scale Density Fluctuations in the Solar Wind 
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Abstract. We motivate the importance of studying Idnetic scale turbulence for understanding the macroscopic properties of 
the heliosphere, such as the heating of the solar wind. We then discuss the technique by which kinetic scale density fluctuations 
can be measured using the spacecraft potential, including a calculation of the timescale for the spacecraft potential to react to 
the density changes. Finally, we compare the shape of the density spectrum at ion scales to theoretical predictions based on a 
cascade model for kinetic turbulence. We conclude that the shape of the spectrum, including the ion scale flattening, can be 
captured by the sum of passive density fluctuations at large scales and kinetic Alfven wave turbulence at small scales. 
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1. INTRODUCTION 

The solar wind contains fluctuations at a broad range of 
scales: from large scale solar cycle variations down to 
small scale turbulence at plasma kinetic scales. While 
studying kinetic plasma turbulence is of intrinsic interest, 
it is at these scales where plasma heating is thought to 
occur, so determining the nature of this turbulence is also 
important for understanding the macroscopic properties 
of the heliosphere. 

For example, it is well known that the fast solar wind 
proton temperature does not vary with radial distance 
R as expected for isotropic adiabatic expansion. Fig. 1 
shows the radial variation of proton temperature mea- 
sured by Helios Plasma Experiment [1]. The plot con- 
tains only data with a small collisional age Ac < 0.01, 
where Ac is the ratio of solar wind transit time to pro- 
ton collision time, e.g., [2]. This selects the fast, hot, low 
density wind, i.e., the purest examples of collisionless 
"fast wind." The radial power law is -0.69 ± 0. 17, which 
is significantly shallower than for isotropic adiabatic ex- 
pansion, for which the power law is -All> for an adia- 
batic index of 7= 5/3 [3]. This non-adiabatic expansion 
is well known both inside [4-6] and outside [6-9] 1 AU. 
In a collisionless plasma, one should ideally consider the 
parallel and perpendicular temperatures separately, and 
they also display non-adiabatic behavior [4, 10, 11]. 

Dissipation of plasma turbulence is a prime candidate 
for the additional heating required for the non-adiabatic 
radial temperature profiles in the inner heliosphere [12] 
(in the outer heliosphere, heating from pickup ion gener- 
ated waves is thought to dominate [13]). In order to un- 
derstand solar wind heating, therefore, we need to know 
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FIGURE 1. Variation of proton temperature with heliocen- 
tric distance for collisionally young (A^ < 0.01) solar wind. The 
darkness represents the number of points in each bin, normal- 
ized to the maximum number for each R. The best fit line (red) 
to the peaks of a Gaussian fit at each R is significantly shallower 
than for isotropic adiabatic expansion (blue line). 



how the turbulence dissipates at kinetic scales. While 
there has been an increasing number of measurements of 
turbulence in this range, there is currently still disagree- 
ment about its nature (see, e.g., [14-16] and references 
therein). 

Recently, Chen et al. [16] measured the density fluc- 
tuation spectrum of solar wind turbulence between the 
ion and electron kinetic scales, finding a spectral index 
of —2.75 ±0.06. Here, we investigate this topic further 
by examining the nature of the flattening of the density 
spectrum at the ion scales. 



2. SPACECRAFT CHARGING 



The density spectrum in Chen et al. [16] was measured 
using the spacecraft potential [17] of ARTEMIS [18] 
as a proxy for density [19]. This technique requires the 
spacecraft potential to react quickly enough in response 
to plasma density changes. In this section, we describe 
the spacecraft charging, including a calculation of the 
charging timescale. 

The charge on a spacecraft is given by Q = CVsc, 
where C is the spacecraft capacitance and Vsc is the 
potential of the spacecraft with respect to the plasma. The 
time dependence, therefore, is given by 



dV,c/dt = /t/C, 



(1) 



where It = dQ/dt is the total current to the spacecraft. In 
typical sunlit conditions, the dominant contributions to /t 
are the flow of electrons from the plasma to the space- 
craft from thermal motions /pi = —n^eA^k^T^/ {Inme) 
[20] (rig = electron density, e = magnitude of electron 
charge, A = spacecraft surface area, = electron temper- 
ature, me = electron mass), and the flow of photoelectrons 
from the spacecraft to the plasma Ipee^^""^^'"', where /pe 
is the photoelectron current at Vsc = and Tpe is the typ- 
ical photoelectron energy in eV. We ignore higher order 
effects, e.g., focusing of the thermal electrons, ion cur- 
rents, probe bias currents, etc. The total current is 



It 



(2) 



We now examine the spacecraft potential reaction to 
density fluctuations. Separating the time varying quanti- 
ties into mean and fluctuating parts, /pi — /pio + 5/pi and 
Vsc = Vsco + 5Vsc, and inserting Eq. 2 into Eq. 1 we obtain 



ddVs, 
dt 



Ipio + 5lpi + Ip,e-^''^^»+^''-^/^^^ 



(3) 



Subtracting the equilibrium equation, in which the cur- 
rents are balanced (0 = /pio +^pee^^'^'"^^P'=), and assuming 
5Vsc ^ Tpe, Eq. 3 becomes 



dSVs, 
dt 



5/pi + ^5y. 



(4) 



Eq. 4 is a first order linear differential equation, which 
can be solved to find the time dependence of 5Vsc given 
an instantaneous change in plasma current 5/pi caused 
by a change in the plasma density 5«e (note that in the 
solar wind, temperature fluctuations can be neglected 
[19, 21]). The solution to the equation is 
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(5) 



where T = — Crpe//pio (note that /pio is negative so T 
is positive). This describes exponential relaxation of 
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FIGURE 2. Magnitude and phase of frequency dependent re- 
sponse of spacecraft potential fluctuations (5 Vsc) to fluctuations 
of thermal electron current to the spacecraft (dipi). 



the spacecraft potential to the new equilibrium Vsco + 
7pe5/pi//pio with time constant T. An increase in elec- 
trons flowing to the spacecraft will result in an equilib- 
rium with a smaller spacecraft potential. 

We can now make an order of magnitude estimate 
of T for ARTEMIS in the solar wind. Approximating 
the spacecraft by a conducting sphere of radius L « 0.6 
m (which gives the same surface area as the spacecraft 
dimensions 0.8 x 0.8 x 1 m [22]) gives a capacitance 
C = An£\)L « 66 pF. The characteristic photoelectron 
energy is Tpe ~ 1.5 eV [23]. The surface area is A « 
4.5 m^, «e ~ 10 cm^^ and « 10 eV, giving /pio « - 
3.8 /X A. These parameters give a charging time of t « 26 
/is, corresponding to a frequency /c = 1/ {2nT) w 6 kHz. 
Since C °^ L, /pio °^ A °^ I? and Tpe is independent of 
L, the charging time depends on spacecraft radius as 
T oc \/L, so larger spacecraft will charge quicker due 
to their greater surface area to collect extra charge. A 
plasma with higher density and temperature, having a 
larger return current, will also lead to faster charging. 

Alternatively, the frequency dependence of the space- 
craft potential response can be determined by substitut- 
ing 8Vsc = dYsce^'^'f and 5/pi = Slpxe^'^'f into Eq. 4: 



(6) 



The magnitude and phase of this response curve is plot- 
ted in Fig. 2. For frequencies / ^ /c, the potential fluc- 
tuations are linearly proportional to and in phase with 
the current fluctuations and do not depend on frequency, 
therefore the density fluctuation spectrum can be well 
measured by a simple calibration curve, as in [16]. For 
frequencies />/c, the response is frequency dependent 
and while the density fluctuations can still be inferred 
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FIGURE 3. Density fluctuation spectra of 17 solar wind 
intervals normalized in scale to the proton gyroradius p;. 



(as long as the potential fluctuations are measurable), a 
correction for the response curve would be required. 

Finally, we note that this derivation requires 8Vsc ^ 
Jpe- For the frequencies / > 10^^ Hz considered here, 
SVsc < 0.1 V so the approximation is well satisfied; 
for larger amplitudes, e.g., at shock crossings [24], the 
response will differ To conclude this section, the density 
fluctuation spectrum of solar wind turbulence for / <C 6 
kHz can be well measured using the spacecraft potential. 



3. ION SCALE FLATTENING 

Fig. 3 shows all 17 of the density spectra discussed in 
Chen et al. [16], with frequencies converted to wavenum- 
ber k under Taylor's hypothesis and normalized to the av- 
erage proton gyroradius of each interval pi. Frequencies 
greater than 15 Hz and less than 5 times the inverse in- 
terval length have been excluded for reliability. The mea- 
sured spectral index of -2.75 for 3 < kpi < 15 [16] is 
marked, as well as a -5/3 inertial range spectral index 
[25]. A slight flattening of the spectrum can be seen in 
between, marked with a -1.1 slope to guide the eye. 

The flattening at 0. 1 < kpi < 1 has been observed pre- 
viously, both at 1 AU [26-29] and in the near-Sun so- 
lar wind [30, 31] and has been attributed to either pres- 
sure anisotropy instabilities [32] or the increased com- 
pressive nature of kinetic Alfven wave turbulence at the 
ion gyroscale [33-36]. In particular, Harmon and Coles 
[35] and Chandran et al. [36] modeled the spectrum as a 
sum of density fluctuations passive to the Alfvenic tur- 
bulence, which dominate at large scales, and active den- 
sity fluctuations from kinetic Alfven wave (KAW) turbu- 
lence, which dominate at small scales. At the crossover 
point, a flattening is naturally obtained with a shape that 
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FIGURE 4. Density spectra (solid blue) at ion kinetic scales 
for (a) high jSj and (b) low jSj. Theoretical spectral shapes 
(dashed black) with their passive (dotted green) and active 
(dotted orange) components are also shown. 



depends on plasma parameters. 

Fig. 4 shows the density spectra for the two intervals 
with the highest and lowest values of proton beta, j3i = 
3.61 and 0.31. The proton to electron temperature ratios 
are similar for these intervals, Ti/Te = 0.74 and 0.76, and 
are typical for slow solar wind. Theoretical curves, con- 
structed using the technique of Chandran et al. [36] for 
the measured parameters, are also marked. The curves 
come from a kinetic turbulence cascade model [37] and 
consist of a passive contribution, which scales like the 
perpendicular magnetic field and an active contribution, 
calculated from the KAW eigenfunctions. To determine 
the relative amplitudes of these contributions, the param- 
eter F = [(5n/«o)/(5vx/vA)]^ (see Eq. 1 of [36]) was 
measured for each interval from the density and veloc- 
ity spectra at frequencies 10^^ Hz < / < 5 x 10^^ Hz. 
Note that the total power normalization of the theoreti- 
cal curves is arbitrary; they have been plotted above the 
measured curves for clarity. 

It can be seen from Fig. 4 that the shapes of the the- 
oretical curves are in close agreement with the measure- 
ments, except, perhaps, for the slope at high frequencies. 
In particular, the flattening at ion scales is well captured, 
with the lower j3i interval showing a more prominent flat- 
tening. This is consistent with the increased compressive 
nature of KAWs at low j3i and is also consistent with 



the very prominent flattening seen in near-Sun measure- 
ments, where j3i is very low [30, 31]. This effect may also 
explain the density spectra of Kellogg and Horbury [29], 
which show more prominent ion scale flattening at low 
density, i.e., periods likely to have a low j3i. 

It is also interesting to note how F varies with j3i. Fig. 4 
shows that the lower J3i interval has a larger F, meaning 
that the relative passive density fluctuations are larger. 
There are two possible explanations for this. Firstly, it 
is thought that the passive density spectrum above ion 
scales consists of kinetic slow mode like fluctuations 
[38, 39], in which the density fluctuations are larger at 
lower j3i. Secondly, the compressive fluctuations are less 
strongly damped at lower j3i [40, 41]. 

Finally, we note that other explanations for the flatten- 
ing [31, 32, 35] and kinetic turbulence in general have 
been suggested. While there is no space here to discuss 
these possibiUties, the model that we consider [36] pro- 
vides a good match to the current observations. Since the 
flattening is expected to be more prominent for lower j3i, 
it should become more easily detectable in situ with fu- 
ture missions Solar Orbiter and Solar Probe Plus as they 
travel in closer to the Sun. 
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